Neutron EDM in Four Generation Standard Model 
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A fourth generation of quarks, if it exists, may provide sufficient CP violation for the baryon 
asymmetry of the Universe. We estimate the neutron electric dipole moment in the presence of a 
fourth generation, and find it would be dominated by the strange quark chromoelectric dipole mo- 
ment, assuming it does not get wiped out by a Peccei-Quinn symmetry. Both the three electroweak 
loop and the two-loop electroweak/one-loop gluonic contributions are considered. With mv, rn t > at 
500 GeV or so that can be covered at the LHC, and with a Jarlskog CPV factor that is consistent 
with hints of New Physics in b — > s transitions, the neutron EDM is found around 10~ 1 e cm, still 
far below the 10 -28 ecm reach of the new experiments being planned or under construction. 

PACS numbers: 13.40.Em 11.30. Er 14.65. Jk 



I. INTRODUCTION AND MOTIVATION 

The Kobayashi-Maskawa (KM) mechanism [l| for CP 
violation (CPV) has been verified by the Belle and BaBar 
experiments [2] . Constituting the flavor and CPV part of 
the Standard Model (SM), it falls short of the observed 
Baryon Asymmetry of the Universe (BAU) by many or- 
ders of magnitude. However, it was pointed out that, by 
extending to four quark generations [3|, SM4, the KM 
picture may have enough [4j CPV phase for BAU. The 
strength of phase transition, needed to satisfy the third 
Sakharov condition, i.e. departure from equilibrium, re- 
mains an issue. But interest has renewed [5j in the direct 
search of fourth generation quarks at hadron colliders, 
where the LHC should finally be able to discover, or rule 
out once and for all [6], fourth generation quarks. 

The long quest for neutron electric dipole moment 
(nEDM) has been motivated by BAU, as the latter im- 
plies the existence of new CPV sources beyond SM. Given 
the large jump in CPV, it is of interest to ask what nEDM 
value one might expect for SM4. The KM mechanism 
cleverly shields itself from nEDM. At the one weak loop 
level, the CKM factors come always conjugate to each 
other so the CPV phase cancels. It was shown [7| by Sha- 
balin, at two loop in SM, that the sum over all diagrams 
for single quark electric dipole moments (qEDM) van- 
ish. It was then shown that, bringing in a further gluon 
loop (two-loop electroweak/one-loop gluonic) breaks the 
identical cancelation, leading to dd ~ 10~ 34 ecm [8|, [9j. 
Considerations of long distance (LD) effects allow a value 
of d n that is two orders of magnitude higher [lCf . 

The current limit for nEDM, 2.9 x 1(T 26 ecm at 90% 
C.L. |ll| . is from the RAL-Sussex-ILL experiment which 
operated at Grenoble. Compared with dropping an or- 
der of magnitude per decade [12J since the 1950s, the 
pace has slowed. The chief limiting factor is the num- 
ber of neutrons. There is, however, a renewed effort, 
by the CryoEDM collaboration at Grenoble, the nEDM 
collaboration at PSI, and the nEDM collaboration at 
the Spallation Neutron Source (SNS) at Oak Ridge Na- 



tional Lab, to push [12, [T^ nrst towards O(10 27 ) ecm, 
then eventually down to 10 -28 ecm. There are efforts 
also at J-PARC [Hj], and TRIUMF [H[, with similar 
aims. All these efforts use ultra cold neutrons (UCN), 
hence many more usable neutrons, together with stronger 
electric fields, and better magnetic field control. The 
SNS experiment is the most innovative, injecting polar- 
ized 3 He into liquid 4 He. The 3 He serves as both [12[ 
"co-magnetometer" and as detector for neutron capture, 
while the liquid 4 He serves as scintillator, as well as a su- 
perthermal source for UCN, and it is much more tolerant 
of a higher electric field. 

These new ambitious experiments are motivated in 
part by the cosmological limit of 10 -28 e cm, if super- 
symmetry (SUSY) is relevant [ll| for BAU. What is the 
situation for SM4, the Standard Model with four quark 
generations, given that CPV is greatly enhanced? In this 
paper we address this issue. 

The nEDM for SM4 was already considered in the 
1990s (l7Hig |. However, they went largely unnoticed be- 
cause the fourth generation fell out of favor in a similar 
time frame, by the twin issues of neutrino counting and 
electroweak precision tests [2, 3]. But since the discov- 
ery of atmospheric neutrinos, we know the neutrinos have 
mass, which calls for New Physics beyond SM. The fourth 
generation neutral lepton [3J, which does not enter our 
discussion, if it exists, must be heavy. More recently, it 
was pointed out [20j that the electroweak S and T vari- 
ables do allow a somewhat split (less than Mw) but close 
to degenerate fourth generation, if one allows the Higgs 
boson to be heavier at the same time. 

In the following sections, we first collect some relevant 
formulas, then proceed to make a numeric estimate. We 
end with some discussions and a conclusion. 



II. SOME RELEVANT FORMULAS 

Starting from the effective Lagrangian of all CPV op- 
erators up to dimension 5, the neutron EDM was evalu- 



ated in the QCD sum rule framework. In terms of quark 
EDMs and chromo-EDMs (CEDM), one has [n| 
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where 1/m* = l/m„ + l/m d + l/m s = l/m u + l/m d , 
9 = J2q@q + ®G is the combined quark and gluonic 9 
term, a^ — e q (2n ± £), and x, m§, k and £ are conden- 
sate susceptibilities. The large factor of 3 uncertainty in- 
herent in the overall 0.4 ±0.2 coefficient reflects the large 
hadronic uncertainty, as determined in the sum rule ap- 
proach. Thus, our estimates that follow are only aimed 
at the order of magnitude. 

The interesting subtlety is that, when a Peccei-Quinn 
symmetry [22| is invoked to remove the 9 term (setting it 
to zero), it induces additional CPV terms [23| to the ax- 
ion potential that completely cancels the strange quark 
CEDM (sCEDM) contribution [Hj], independent of de- 
tails of the axion potential. While remarkable, as we 
shall see, the sCEDM is of the greatest interest in SM4. 
Furthermore, three decades of axion search has so far 
come to naught. Given that there are models of spon- 
taneous CPV, such as the Nelson-Barr mechanism [25| . 
where 9 = while the KM phase (phases for SM4 case) is 
generated by spontaneous symmetry breaking, we shall 
ignore the 9 term while keeping the qCEDM terms. 

We follow Ref. [2l| and take the numerical values m 2 , = 
0.8 GeV 2 , x = -5.7 ± 0.6 GeV -2 , £ = -0.74 ± 0.2, and 
K = —0.34 ± 0.1. Numerically, we then have 

d n = (0.4 ±0.2) 1.9 x lO- 16 0ecm-O.O8ed s 

+ 1.8e d d - 1.4e d u + {Ad d - d u ) . (2) 

The strange quark CEDM entered by m s being light 
enough, such that it partakes in the qq condensation. 

Analyzing the flavor structure of a typical three loop 
diagram shows why the strange CEDM is highlighted, de- 
spite a smaller coefficient in Eq. ^. A typical three loop 
diagram involves two nonoverlapping W boson loops, 
with one Z or gluon loop. Following the quark line, the 
/ quark (C)EDM has the following form 0, [H| 
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j,k,l 

= l -J2 Im (VffVjkVttVtf) f(jk I - lkj)f, (3) 

3,k,l 

where f,j,k,l stand for both flavor indices and the cor- 
responding Green function. The antisymmetry in Eq. ([3]) 
is at the root of Shabalin's argument. 

Since we shall consider rather heavy t' and b' quarks in 
the loop, typical loop momenta would be at these large 



values. Therefore, one can take c = u=u, d = s = b = d 
as all effectively massless in loop propagators. One then 
easily sees that the (C)EDM of the u quark vanishes. 
That is, performing the sum over j and I in Eq. (|3|) for 
/ = u, using the unitarity relation V* d V kd + V* s V ks + 
V* b Vkb = S u k - V* b ,Vkb' and the "degeneracy" in mass 
for the d, s, b propagators, one gets 



i^Mv^VkjVjWJujkh 



j,k,l 

= ^ lm ( V ub'Vkb>V k \,V ub ,)u(dkb' - b'kd)u = 0, (4) 

k 

as phases pairwise cancel. Effectively there are only two 
generations in the loop. 

The case for / = d, s is therefore more interesting. By 
similar token, one has 

j,k,l 

= iIm(V t * f V tb V t 1 b V t , f )f[t (d - b') t' - t'(d - b') t 
+ t'(d-b')u-u(d-b')t' 
+ u{d-b')t-t(d-b')u]f. (5) 

In Eq. ([S]), we have spelled out the sum over k, after 
utilizing CKM unitarity as before. The sum of Green 
function factors contain the "degeneracy" of c = u = u 
and d = s = b = dm the loop. But one should treat the 
CKM coefficient with care, where use has been made of 
the rather good approximate relation [17| 

lm(V t }V tb ,V t 1 b ,V t , f ) - -lm(V t }V tb V t r b V t/f ) = J f , (6) 

which is a consequence of the smallness of V ub (assum- 
ing that other CKM elements that enter are not much 
smaller). As noted in Ref. (27|, for / = s, this relation 
effectively means that the CKM "triangle" (degenerated 
from a quadrangle) governing b' —¥ s transitions have es- 
sentially the same area as the b — > s triangle. We will 
turn to numerical analysis in the next section, but we can 
already see that the CKM factor for / = s is much larger 
than for / = d, which is the reason why we highlight the 
strange quark CEDM. 

The s quark CEDM arising from the two-W loop plus 
one gluon loop diagram was estimated [17J j using the ex- 
ternal field method, to double log accuracy and in the 
large Nc limit, with the result of 
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Replacing the gluon by a Z boson loop, one might 
think it should be subdominant. But since there are 
rather large quark masses in the loop, it implies that 
rather large loop momenta might be relevant. The 
derivative coupling nature of the longitudinal Z boson, 
or equivalently, that the Goldstone boson couples to the 
heavy quark masses hence is nondecoupled, effectively 
voids the above intuition. By an ingenious argument of 



limiting to large loop momenta and involving longitudi- 
nal vector bosons, the authors of Ref. [18j were able to re- 
duce the three-loop calculation effectively to calculating 
three one-loop integrals, and the core of it is an effective 
i — > fZ transition involving the heavy fourth generation 
quark in the loop. This is the familiar Z penguin [28l.l29l]. 
and indeed it has been found [6] that b' — > bZ and b' — s- bg 
transitions are not too different in strength. The upshot 
of the estimate (with the brutality of setting all loga- 
rithms to order 1) of Ref. [18|], done in f996, is 



d{ z) = -J s m. 



G 



F Oi w 



V2 (4tt)' 



4M 4 , 



i I m t' 



(8) 



Comparing Eqs. ([7]) and ([8]), aside from the double 
versus single logarithm, one can see from aw/Myy = 
\J~2Gf /tt = 1/nv 2 that one is comparing 5N c a s /6 with 
Aj//47r. The gluonic effect is enhanced by the color factor, 
but the Yukawa coupling grows with m 2 / . Compared lit- 
erally, they are actually comparable. On the other hand, 
in arriving at Eq. ([8]), one has set all logarithms to f. 
In this spirit, both the double log (including the 1/2!) in 
Eq. (|7|) and the single log in Eq. (j8|) should be treated 
as order one. Then, the gluonic effect would be subdom- 
inant to the Z effect, for t 1 and b 1 masses of order 500 
GeV (or higher), a nominal value used by Ref. [18(, and 
which we shall use in the next section. 

Given the roughness of these calculations, and the 
great difficulty in calculating genuine three electroweak 
loop diagrams, we shall take the estimate of Eq. © for 
our subsequent numerics. 



III. NUMERICAL ESTIMATE 

We shall use m t < ~ m b > ~ 500 GeV as our nominal 
fourth generation quark mass. Above this value, one 
would pass through the unitarity bound [30], and the 
Yukawa coupling starts to turn nonperturbative. For the 
other quark masses, we take m u = 2.5 MeV, rrid = 5 
MeV, m s = 100 MeV, m c = 1.3 GeV, m b = 4.2 GeV, 
and mt — 165.5 GeV, where heavy quark masses are in 
MS scheme. The light quark mass values were in fact 
implicit in Eq. ©. For the CKM products J s and Jd-, 
we take the nominal fit [3lJ to flavor data performed for 
mt' ^ 500 GeV, where V t , b ~ -0.1, V Vs ~ -0.06 e - l75 °, 



and Vf 
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-0.003 e 
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, we get 

Im(V t * s V tb V t r b V t , s ) 

Im(V t * d V tb V t 1 b V t , d ) 



2.4 x 10~ 4 , 
: 1.7 x 10~ 7 . 



(9) 
(10) 



Note that J s could be measured [32J in the next two years 
at the LHC, but Jd would be harder to disentangle. 

Although the range of uncertainty is large, it is clear 
that | J s | ^> \Jd\, which correlates with the fact that 
b — > d transitions, including in Bd mixing, show little 
sign of deviation from SM expectations, while for b — s- s 
transitions, and especially in B s mixing, we have several 
indications for sizable deviations from three generation 



SM. Note that the study of Ref. [31] predated the summer 
Tevatron update on sin2$s s , the CPV phase of the B s 
to B s mixing amplitude, and predicted a value lower than 
the -0.5 to -0.7 value given [H| in 2007 for m v ^ 300 
GeV. In any case, for our purpose, we see from Eqs. (|9]) 
and (ITUl) that, despite the smaller coefficient in Eq. ©, 
the sCEDM is the most important quark contribution to 
d n , because of the CKM factor. 

Putting in numbers, we find from Eq. ([8]) that 

J< 4) ~-4x 10~ 16 GeV" 1 ~ -0.8 x 10~ 29 cm, (11) 

where the W-W-g 3-loop effect of Eq. ([7]) is treated as 
subdominant. Treating the sCEDM as the leading effect 
in Eq. ©, then 



4 4) = (2.2 ± 1.1) x 10" 31 ecm, 



(12) 



where the superscript indicates that this is the estimated 
effect from the fourth generation. 



IV. DISCUSSION AND CONCLUSION 

Though the estimate of Eq. (fT2")l is considerably larger 
than the SM result with three quark generations, be it the 
quark level [9j , or hadronically enhanced [10( , it is very far 
from the sensitivities of next generation experiments [l2|- 
Ua ]. nominally at 10~ 28 ecm. What is worse, or perhaps 
intriguing, is that the sCEDM effect of Eq. (Ti~Tj) is at the 
mercy of the Peccei-Quinn symmetry. If a PQ symmetry 
is operative in Nature, then the sCEDM effect is precisely 
canceled [23, [24|. In this case, one has a reduced formula. 
Rather than Eq. ©, one gets [2l| 

rfP Q = (0.4 ± 0.2) [l.6e(2(& + d u ) + (4d d - d u )] ,(13) 

i.e. only dependent on the naive constituents of the neu- 
tron. Note that the axion potential has also modified the 
d and u quark CEDMs, as it brings in analogous terms 
to the one that canceled away the sCEDM effect. 

From Eq. (|4|), d u and d u vanish. Thus, we need only 
consider dd and dd for SM4. Assuming again that the 
analogue of Eq. ((5} dominates over the gluonic counter- 
part, we use the formulas of Ref. [21( for d s and Eq. ([8]) 
to obtain dd and dd by simply shifting the CKM index, 
i.e. shifting from J s of Eq. © to J d of Eq. (fTUD, and 
replacing m s by irid- We find 

7 (4) 



-3xl0~ 34 cm, d d 4) ~ -4x 10~ 34 ecm, (14) 



d (4) 
and, from Eq. (|T3| . 

4 4 ) p Q = -(l±0.5)xlO- 33 



ecm, 



(15) 



where dd contributes roughly twice as dd- These should 
be taken as very rough estimates. Note that, unlike 
Eq. ©, it would take some time to refine Eq. (fT0"|) . 

We see that d d is stronger than in SM Q, while 
d n , as estimated through sum rules, is at 10~ 33 ecm 



level. Other LD effects might bring about another [10| 
order (or maybe two) of magnitude enhancement. But 
it should be clear that, with PQ symmetry operative, 
even with four quark generations, nEDM is way below 
the sensitivity of the next generation of experiments. 

But what about LD effects to the sCEDM-driven result 
of Eq. (fTi-j) . without PQ symmetry? Operators beyond 
dimension 5, such as dimension 6 operators involving two 
quarks in the neutron, are beyond the scope of our in- 
vestigation. Even with WWg loops, Eq. (|7|), of similar 
order of magnitude and constructive with the WWZ loop 
result of Eq. flSJ} (hence Eq. ([IT])), df 1 without PQ sym- 
metry still seems a couple of orders below 10 -28 e cm. 

One could also estimate from s to d transition elements 
induced by gluon or Z , enhanced by p ion-nucleon cou- 
pling at long distance. It is possible 34[ that Kl — > ifivv 
gets enhanced by a factor of 100 in SM4, which means 
a factor of 10 at amplitude level. Given the two orders 
of magnitude LD enhancement in SM [lOj, it could be 
brought up by another order of magnitude in SM4, hence 
to 10~ 31 e cm. This is rather similar to the result in 
Eq. (|12[) . We remark that, one way to probe the sEDM 
and sCEDM effects, even if PQ symmetry is operative, 
would be to measure hyperon EDM, which is very inter- 
esting in its own right. A rough estimate, by extrapolat- 
ing from Eq. (fT5)l. gives the order at 10 _29 -10 -28 ecm. 
But if 10 _28 -10 -27 e cm is still the challenge of the next 
decade for neutron EDM experiments, the measurement 
of the EDM of the less abundant, shorter-lived hyperons 
would be still farther away. 



The CMS experiment has recently stated that "data 
exclude SM4 Higgs boson with mass between 120 and 600 
GeV at 95% C.L.", while the t' quark is not seen below 
450 GeV [35]. One may be entering the strong coupling, 
beyond the unitarity bound regime of m t ' > 500 GeV, 
with an associated composite, massive (hence broad) 
"Higgs" boson. In this regime, the short distance results 
of Eqs. (J7J and © start to fail even as rough estimates, 
as the loop functions themselves turn nonperturbative. 

In conclusion, with four quark generations and with 
Peccei-Quinn symmetry operative, the neutron EDM is 
slightly enhanced above the SM value, but no more than 
an order of magnitude, hence much below the sensitivities 
of the next generation of experiments. If PQ symmetry is 
absent, then a large enhancement is possible through the 
s quark CEDM, which is correlated with possible effects 
in b —¥ s transitions that are of current interest. However, 
it is still unlikely that the neutron EDM could reach the 
10~ 28 e cm level sensitivity that may be probed during 
the next decade. 
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